Universality of Anderson transition in two-dimensional systems of symplectic 

symmetry class 
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We investigate localization of noninteracting particles with spins higher than 4 in a two- 
dimensional random potential in presence of spin-orbit coupling. We consider an integer spin (s = 1) 
and a half-integer spin (s = §) belonging to orthogonal and symplectic symmetry classes, respec- 
tively. We show that particles with integer spin are localized and those with half-integer spin exhibit 
Anderson transition. The transition belongs to universality class of conventional symplectic model 
for spin-| particles. 
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I. INTRODUCTION 

Symplectic class has a rich physical content among 
the symmetry classes in the Wigner-Dyson classification 
of random matrices. Spin-orbit scattering provides the 
common physical realization of this class. From the sym- 
metry point of view such system is invariant under time 
reversal (7~) but not under spin rotation (S). For a sys- 
tem with non-integer spin we have T 2 = — 1 and the wave 
function (WF) has a rotational periodicity of 4ir therefore 
in average, time-reversed paths in the multiple-scattering 
picture interfere destructively. 1 As a result of the de- 
structive interference, conductivity is enhanced and we 
have weak anti-localization rather than weak localization. 
This is one of the mechanisms of criticality in two spatial 
dimensions which is specific for symplectic class in the 
framework of early Wigner-Dyson classification. 

Several other realizations of symplectic class have been 
identified which exhibit distinct universal behavior. For 
this class the homotopy group of a model manifold is 
nontrivial so the a model action allows for inclusion of a 
topological 9 term (see Ref. 2 for a recent review). Such 
topological term is responsible for quantum Hall critical- 
ity in the unitary class. Spin-orbit coupling have pro- 
vided similar topological phases in the T-invariant sys- 
tems which has been the subject of an intense activity 
in recent years. 3 Remarkably, quantum spin Hall (QSH) 
which is a novel phase induced by Z2 topological term 
in the symplectic ensemble. A new universality class 
of Anderson transition emerges in presence of this topo- 
logical structure between the metallic and QSH phases. 4 
It should be mentioned however that the corresponding 
Chalkcr-Coddington network model which allows to have 
odd number of Kramers doublets (i.e., nontrivial topol- 
ogy) does not capture this critical behavior. 5 Another 
realization of symplectic class with nontrivial topology 
appears in a two-dimensional system of Dirac fcrmions 6 ' 7 
which yields a unconventional scaling f3 function. Two 
different scenarios are proposed one of which predicts an 
extra attractive fixed point in the strong-coupling limit 6 
and the other one implies on derealization of all states 
even in strong disorder limit. 7 Based on semiclassical ar- 
guments it is also shown that spin-orbit scattering may 



induce a novel universality class in the regime of integer 
quantum-Hall effect. 8 Whereas in the unitary ensemble, 
presence or absence of spin-rotational invariance does not 
change symmetry class. These examples imply the fact 
that in spite of a complete mathematical classification of 
symmetry classes, universality classes are not recognized 
so far. 

Appreciate to new advances in designing periodic po- 
tentials by standing waves of light, many experiments 
which are not possible to arrange for electrons can be 
simulated with quantum motion of cold atoms. Espe- 
cially it just recently became possible to do careful ex- 
periments on localization of noninteracting matter waves. 
Even more fascinating games could be done by changing 
the polarization of the beams. So the internal degrees of 
freedom of the atom can be coupled to the momentum 
of the beam and produce an effective spin-orbit like term 
in the Hamiltonian. Once this could be done, one can 
search for new universalities and topological properties 
of WFs with higher tunability. 9 

A natural generalization in this direction is to con- 
sider atoms with higher number of internal degrees of 
freedom or particles with higher spins. Here we want 
to address whether higher spins in spin-orbit interaction 
can change the universality class of transition or not. We 
use the transfer-matrix method to calculate the localiza- 
tion length and then extract the critical exponents from 
finite-size scaling analysis. On the other hand multifrac- 
tal spectrum of critical WFs are of universal properties 
hence useful to describe the transition. We will examine 
multifractal properties of higher spin model in compari- 
son with spin-i case. 



II. MODEL 

Several models have been proposed to study An- 
derson localization problem in presence of spin-orbit 
scattering. 10 Regardless of microscopic details, they 
present the same universal features. Spin-relaxation 
length is an important irrelevant length scale in these sys- 
tems. SU(2) model 11 has smaller spin-relaxation length 
since spin rotation operators in each link of its lattice, 
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are uniformly distributed. As a result it needs small cor- 
rection to scaling. Having small size effects this model 
provides more accurate calculation of critical exponents. 
Here we use the generalization of this model to describe 
particles with higher spins. We start with the following 
Hamiltonian which is proposed for spin-i particles. 

Here we will let the hopping matrices act on spinors 
of higher rank. So spin indexes (cr, a') take the values 
(s, s — 1, • • • , — s) for particles with spin s. Latin indices 
denote nearest-neighbor sites on square lattice. Random 
on-site potential is distributed uniformly in the interval 
[ — 'T' If]- Energy scale will be set by V = 1. 

TZ iJ, s are (2s + l)-dimensional irreducible repre- 
sentation of SU(2) group. In terms of Euler an- 
gles they have the following description TZ^ , = 
D%A<*ii>Pii>Hi) = e i (-«+^) d « (T ,( A .) ) whcrc 
d s aa ,{Pij) = {stj\er tl3iiS v/ h \sa') is the matrix element of 
rotation operator around y axis. Angles a, f3, and 7 are 
distributed randomly in different links of lattice such that 
rotation matrices IZ 13 have uniform distribution with re- 
spect to the Haar measure on SU(2) group. Namely, a 
and 7 are distributed uniformly in the interval [0, 2n) 
and /? is chosen from interval [0, ^] with distribution 
PCS) = sin(2/3). 

III. LEVEL STATISTICS 

Regarding broken spin-rotational symmetry of Hamil- 
tonian (1), for integer spins it does belong to the orthog- 
onal ensemble. Only for half-integer spins it falls into the 
symplectic ensemble. This is known from Wigncr-Dyson 
classification of random matrices. 12 An essential differ- 
ence in the spectrum of two cases is Kramers degeneracy 
of energy levels of later which is robust against disorder. 

In this section we demonstrate above-mentioned rela- 
tion with the Wigner-Dyson symmetry classes numeri- 
cally. The simplest quantity which can be used to de- 
termine the statistical properties of energy levels is the 
nearest level-spacing distribution. To have comparable 
results with the distributions of random matrix theory 
we look at the distribution of 5 n = which 
is unfolded level spacing. Denominator is the ensem- 
ble average of level spacing which is proportional to in- 
verse density of states (DOS). So unfolding procedure 
is needed when the DOS has large variations within 
the energy range under consideration. Obviously first- 
order moment of distribution function P(S) is fixed, i.e., 
(S) = J Q 5 P(S)d5 = 1 for unfolded spectrum. We de- 
termine the distribution function P(S) for two s = 1 and 
s = I cases in the metallic regime, after removing the de- 
generacy in the later case. By metallic regime we mean 
weak disorder for which WFs have large overlap and com- 
parable localization length with the system size. We find 




FIG. 1: (Color online) Level-spacing distribution function for 
disorder strength W = 0. Spin l(GOE), spin §(GSE) and 
corresponding Wigne-Dyson distributions with /? = 1,4, re- 
spectively. Inset: log-log plot which shows power-law behav- 
ior at S —¥ 0. Line segments represent functions proportional 
to 5 and S 4 . 

good agreement with GOE (j3 = 1) for s = 1 and GSE 
(/3 = 4) for s = § cases. The results which are shown in 
Fig. 1 are obtained by diagonalizing 10 3 Hamiltonians of 
lattice size 20 2 and disorder width W = 0. For nonzero 
but small values of W also we obtain the same results. In 
the strong disorder limit which all states (for both cases) 
tend to be localized one naturally expect to see Poisson 
distribution. Difference between two cases would reveal 
in the intermediate disorder strength. 

IV. TRANSFER MATRIX 

To find a precise insight in to the localization prop- 
erties of these models and to explain the differences in 
the thermodynamic limit we implement a finite-size scal- 
ing analysis. In the following we will study renormal- 
ized localization length (RLL), A = -j^, on the quasi- 
one-dimensional geometry, where A m is the localization 
length on strip. We utilize transfer-matrix method 13 
to calculate the minimal Lyapunov exponent, inverse 
of which is the largest length scale of spatial exten- 
sion of wave function. Dimension of transfer matrices 
for spin s is TV = 2 (2s + 1)M with M being the sys- 
tem size in the transverse direction. Lyapunov expo- 
nents appear in (—7,7) pairs for integer spin case due 
to symmetry of transfer matrices. Furthermore for half- 
integer spin case each 7 appears twice due to Kramers 
degeneracy. According to these symmetries we need to 
evolve = (2s + 1)M vectors of length TV to calcu- 
late minimum positive 7. Components of vectors are 

V( 2s +i)j+m = il>2' 3 , V(2 S +i)(j+jvf)+m = where 
n is the number of layer (here a chain of length M), 
j = 0, . . . , M — 1 is the coordinate in the transverse di- 
rection and m takes values 1, . . . , (2s + 1) corresponding 
to cr = s,s — 1,..., — s, respectively. Gram- Schmidt or- 
thogonalization is implemented after each four steps. 
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FIG. 2: (Color online) Left: RLL as a function of energy of 
spin-1 particle with zero on-site disorder [W = 0), right: spin 
| and spin | (for single size M = 8) with W = 1. Inset shows 
a zoom in around the crossing point. 



Let us start with spin-1 particle and zero on-site disor- 
der (W = 0). In the left panel of Fig. 2 we observe that 
A decreases by increasing the size M in the whole energy 
range. We can conclude that states are localized even for 
zero on-site disorder. In other words the randomness in 
spin rotation in passing through different links is enough 
to localize the particle. Results for nonzero on-site dis- 
order are the same and we will not present them here. 
This is what we expect for a system in orthogonal (AI) 
symmetry class. That breaking of spin-rotational sym- 
metry in a T-invariant system with integer spin neither 
changes symmetry class nor develops dclocalized states. 
We should comment on the additional symmetry which 
Hamiltonian (1) may have in absence of on-site disorder. 
Using periodic boundary conditions in the transverse di- 
rection and even M the lattice will be a bipartite lattice 
which is shown to have anomalies at zero energy. 14-16 We 
can see in Fig. 2 that at the center of energy band, A re- 
mains almost constant for different sizes which indicates 
a critical state at E — 0. Away from the band center A 
decreases more rapidly by increasing M, which gradually 
leads in creation of a cusp at E = 0. For an odd number 
of channels (M) with free boundary conditions in trans- 
verse direction particularly A — > oo for this zero mode. 
This is proved 15 analytically for coupled one-dimensional 
chains with j3 = 1,2. This critical state will be ruined 
by addition of small on-site disorder which breaks sublat- 
tice symmetry of Hamiltonian. Unlike the spin-1 particle 
the case with spin-| possesses a band of extended states 
for certain values of on-site disorder strength W. As an 
example, results of A for spin-| particle with disorder 
strength W = 1 and sizes M = 8, 16, 24 are shown in 
Fig. 2. The result for spin-i particle is also included for 
comparison at the same on-site disorder strength. Rela- 
tive errors of data arc 1% for M > 8, 0.8% for M = 8, 
and 0.7% for spin-i case. To reach this accuracy the 
length of strips is increased up to 1.3,1.6,2.4,3.6 x 10 6 



FIG. 3: Ratio of RLLs vs W of spin-| and spin-i particles 
for single energy E = 2 and M = 8. 

for spin 1 with M = 8, 16, 24, 32, 1.4, 1.8, 2.8, 3.8 x 10 6 for 
spin | with M = 8, 16, 24, 32, respectively, and 0.9 x 10 6 
for spin-^ case. There can be seen a mobility edge at 
E rj 3.4 where order of symbols is reversed. As usual 
the states in the band edge are localized and the mid- 
band states are extended. The finite band of the bulk 
extended states make the true metallic phase happen in 
the symplectic class. By increasing the disorder strength 
the band of extended states gets narrower and gradually 
collapses at a critical value W c . 

It is worthwhile to compare RLL of spin-| and spin- 
\ cases more closely. The ratio of RLLs vs W of two 
models is plotted in Fig. 3 for single energy E = 2 and 
AI = 8. The energy is chosen away from the band cen- 
ter and edges to ensure the DOS has considerable value 
in both models. In a range of weak disorder strengths 
(W < 4), the ratio is nearly constant and equals 2. This 
is where the hopping term is dominant or comparable 
with on-site term. In strong disorder limit (W 3> 1) the 
hopping term is negligible, therefore spin degrees of free- 
dom would not have considerable effect on localization 
length. Thus one expects the same RLL for both cases. 
That is what which can be seen also in Fig. 3 at large 
W. In the next sections we will discuss critical exponents 
characterizing the universality class of the transition. 

V. SCALING AND CRITICAL INDICES 

Dimensionless quantity A is one of scaling variables 
which is frequently used for numerical analysis. On the 
basis of one-parameter scaling hypothesis, it can be writ- 
ten in the following form: 

where £,(W, E) is the localization length (insulating side) 
or the correlation length (metallic side) of infinite system. 
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It is not the only length scale in this system. We will 
encounter deviations from scaling, Eq. (2), when other 
(irrelevant) length scales are comparable with correlation 
length. Near the mobility edge, £ diverges as £ ~ \E — 
E c \~ v with critical exponent v. The value of A at critical 
point is also a universal constant and independent of W . 

To calculate the critical exponent we take few sets 
of data (A vs E) close to the mobility edge for M — 
8, 16, 24, 32. There are two ways of fitting the data to 
scaling form (2). Since function f(x) is unknown we can 
either use a Taylor expansion and then obtain the coef- 
ficients by fitting, or take one set as f(x) (by interpolat- 
ing), define a suitable residual of curves and minimize it 
by adjusting critical indices. Let the values of A and E 
in the ith set (corresponding to size Mi) be denoted by 
Aij and Eij. 
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FIG. 4: (Color online) RLL of spin- § particle for W 
sizes M = 8, 16, 24, 32. Inset: after rescaling. 
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A possible definition of residual follows 17 
fl 4EEE' > A « - AfcCM 1 /"^ - E c ))\ (3) 



k i^=k j 



where Afc(x) is obtained by transforming the horizontal 

axis of Akj such as E^j — > M]J V (E^j — E c ) and simple 
linear interpolation. The prime on the third sum denotes 
summation over j's which are in the range of definition 
of Afe(x) and N is the total number of such points. For 
each k one set is taken as the reference curve and the 
other curves are supposed to collapse on it by rescaling. 
Function R reaches its minimum value (~ 0.01) at v = 
2.81±0.18and£: c = 3.362±0.014. Error-bars are roughly 
estimated from width of minimum by using approximate 
expressions given in Ref. 17. Reasonable data collapse is 
obtained for these values of critical parameters (inset of 
Fig. 4). 

We tabulate the values of critical exponents in Table I. 
The exponents of spin-g system are also given for com- 
parison. In spite of relatively larger error-bars in the 
spin- 1 case, estimated exponent is close to that of spin-^ 
particle. Our speculation is that two models belong the 
same universality class. 



TABLE I: Critical exponents: results of spin-^ case are taken 
from Ref. 11. 



Model 


V 


A c 


Spin | 


2.73 ±0.02 


1.844 ± 0.001 


Spin | 


2.81 ±0.18 


1.77 ±0.07 



VI. MULTIFRACTAL SPECTRUM 

Fluctuations of critical WFs exhibit universal features. 
So multifractal exponents which describe the distribu- 
tion of WFs can be used to characterize phase transition. 
Moreover multifractality leads to anomalous diffusion of 
wave packets near the mobility edge. This is reflected 
in power-law decay of return probability which is gov- 
erned by one of those multifractal exponents. We want to 
compare the multifractal spectrum of spin-| and spin-i 
particles at the mobility edge. Scaling properties of mul- 
tifractal measure are encoded in the f{a) spectrum. Here 
we use the direct method of Chhabra and Jensen 18,19 to 
calculate this function. By using the box probabilities 
p(l) = J n ^ d 2 r\fy{r)\ 2 , with 0,(1) being a box of linear 
dimension I, and one-parameter families of normalized 
measures Hi(q,l) = Pi(0/Eif>i(0 f° r cacn va htc of q, 
corresponding values of a(q) and f(q) = f[a(q)] can be 
calculated from 



a(q) 



ln(//L) : 



f(q) = 



ln(//L) 



(4) 



In practice, linear fit to numerators vs ln(l/L) can be 
used to calculate a and / as the slope of fitted line. This 
method requires large system sizes to avoid finite-size 
effect. 
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FIG. 5: (Color online) Singularity spectrum of critical WFs of 
spin-^ and spin-| particles with lattice sizes 1500 2 and 800 2 , 
respectively. Line f(a) — ot is tangent to f(a) curve as it 
should be. Horizontal line shows maximum of f(a) which is 
dimension of support (d = 2). The data from Ref. 21 are 
included for comparison. 



These exponents are related to correlation dimension 
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r(g) by a Legendre transform a{q) = dr(q)/dq, f[a(q)] = 
a(q)q — r(q) and generalized dimension D(q) is defined 
through r(q) = (q — l)D(q). 

The singularity spectrum for single critical states of 
spin-i and spin-| particles are shown in Fig. 4. Cor- 
responding lattice sizes are 1500 2 and 800 2 , respectively. 
The eigenstates are obtained via Lanczos algorithm. Pa- 
rameters of former is picked up from phase diagram ob- 
tained in Rcf. 11 {E = 1, W = 5.952) and the spin-| case 
has (E — 3.362, W = 1). Even thought two states are 
distant in parameter space, their singularity spectrum 
are the same within the error bars. Especially we obtain 
(for WFs in Fig. 4) a^ 2 = 2.174 ± 0.005 and D 1 ' 2 {2) = 
1.66 ± 0.04. These exponents and also whole spectrum 
are compatible with previous calculations based on in- 
verse participation ratio analysis. 20,21 For spin-| case 

a 3 Q /2 = 2.164 ± 0.009 and L> 3 / 2 (2) = 1.71 ± 0.05, where 
Q'o = a(0). Error bars are standard deviation of slopes 
in linear fitting (Fig. 5). 

It should be noted that, obtained results depend 
weakly on realization of disorder so ensemble averaging 
will not change them significantly. 



VII. SUMMARY AND FUTURE WORK 

Localization of particles with spin 1 and | in the pres- 
ence of spin-orbit interaction is studied numerically in the 
SU(2) model. In summary, spin-1 particle belongs to or- 
thogonal symmetry class and is always localized. Spin-| 
case exhibits a transition. Implication of finite-size scal- 
ing results and multifractal analysis on critical WFs of 
this model is that the transition belongs to the conven- 
tional universality in symplectic class. 

An interesting direction of future work would be in- 
vestigation of spin dependence of localization length (see 
Fig. 3) within the Dorokhov-Mello-Pereyra-Kumar for- 
malism for disordered wires. 
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